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Abstract 

Multistep dynamical phase transition from the locked to the running state 
of atoms in response to a dc external force is studied by MD simulations 
of the generalized Frenkel-Kontorova model in the underdamped limit. We 
show that the hierarchy of transition recently reported [Braun et al, Phys. 
Rev. Lett. 78, 1295 (1997)] strongly depends on the value of the friction 
constant. A simple phenomenological explanation for the friction dependence 
of the various critical forces separating intermediate regimes is given. 
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I. INTRODUCTION 



The nonlinear response of a system of interacting atoms to a dc driving force has re- 
cently attracted a great interest (see [|T| p~0|| and references therein). The knowledge of the 
microscopic mechanisms for mobility, friction and lubrication processes are very important, 
in particular, for a better understanding of solid friction at macroscopic level, as well as in 
various fields of applied science and technology such as adhesion, contact formation, friction 
wear, lubrication, fracture, etc. 

One generic example represents a layer of atoms adsorbed on a crystalline surface. The 
adsorbate in this case is considered as an atomic sub-system and the remainder is modeled 
as an external potential, a damping constant and a thermal bath. Such a system can 
be treated within the framework of a generalized Frenkel-Kontorova (FK) model JTT|-|T3||. 
When an external dc force is applied to such a system, its response can be very nonlinear 
and complex. By contrast, the driven motion of a single Brownian particle in the external 
periodic potential has been studied in detail and is now well understood ||14|| . If the force 
F is applied to the particle, the total external potential in the direction x of the force is a 
corrugated plane, with a slope F. For small forces the potential has local minima, and the 
particle is locked. The local minima disappear at forces higher then Ffo = Cue /a, where e 
is the amplitude of the periodic potential, a its period, and C a numerical factor depending 
on the shape of the potential. Thus, when the applied force is adiabatically increased, the 
atom passes from the locked to the running state at Ff , and the mobility B = (v)/F 
(where (v) is the drift velocity) reaches its maximal value Bf = l/mn, where m is the 
atomic mass and 77 the friction coefficient. On the other hand, if one decreases the force 
F adiabatically starting from the running state, the critical force F b0 ~ irjy/me/ir for the 
backward transition to the locked state is different owing to inertia of the system. Namely, in 
the underdamped limit (rj ^ uq, where Uo is the frequency of atomic vibration in the external 
potential), the inequality F^q < Ffo holds, and one can observe hysteresis: in the bistable 
region F&o < F < Ff , the particle is either locked or running depending on its initial velocity. 
For a single particle, the bistability disappears at any nonzero temperature. Besides, since 
the "forward" critical force Ff is independent on the friction, and the "backward" force F^q 
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grows linearly with friction, the width of hysteresis vanishes at 77 > rj* = " " 



4a 

The problem of interacting particles in a periodic potential is much more difficult. For 
the overdamped case (77 ^> u ) the nonlinear mobility of the FK model has been studied in 
a number of papers [lHlIi- By contrast, investigations of the underdamped case are very 
limited. In that context, Persson |J observed hysteretic dynamical phase transition, similar 
to the T = one-particle case, in the MD simulation of a 2D system of interacting atoms 
subjected to a periodic potential. Besides, our recent work [|1(J on the underdamped gen- 
eralized FK model revealed strong collective effects in the dynamics of the dc-force driven 
layer of atoms. When the external force increases, the FK system exhibits a complex hier- 
archy of first-order dynamical phase transitions from the completely immobile state to the 
totally running state, passing through several "slip-stick" intermediate stages characterized 
by the running state of collective quasiparticle excitations of the FK model known as kinks 
||10|| . The scenario of these intermediate transitions depends on whether the concentration 
of atoms corresponds to trivial or complex atomic structure. All the observed transitions 
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are hysteretic and, it is remarkable, that, by contrast with the case of noninteracting atoms, 
the hysteresis survives at nonzero temperature of the system. 

However, the results [|10| have been obtained for one value of the friction constant rj only. 
That is why, in the present work we are interested in the question how the observed dynam- 
ical transitions evolve when the friction changed. We calculate dynamical "phase diagrams" 
of the FK system in the (F,rj) plane for two generic atomic concentrations (Sect. We 
show that several critical forces, separating intermediate stages during the transition from 
the locked state to the running state of the atoms are friction dependent, and we propose a 
simple phenomenological approach (Sect. [TVp which allows to explain these dependences. 



II. THE MODEL 

The detailed description of the generalized FK model under study, of the numerical 
procedure, and also the substantiation of the choice of model parameters can be found in 
Ref. ||10|| : here we only outline the main aspects of the model. The atomic motion is governed 



by the Langevin equation, 



d 

mxi + mnxi + - — 



V sub (xi,yi,Zi) + Vint(\ri - rj\) 



F^ + 6F t {x \t), (1) 



for the x coordinate of i-th atom, and similar equations for y and z. Here m is the atomic 
mass, V su b(x,y, z) the external substrate potential, Vi nt (r) the potential of pairwise interac- 
tion between atoms, r] the viscous friction constant, which corresponds to the rate of energy 
exchange with the substrate, F = {F, 0, 0} the dc driving force, and 5F the Gaussian random 
force with correlation function 

(5Ft\t) 5Ff\t')) = 2 V mk B T5 al3 5 tJ 5(t - t>). (2) 

For a better representation of the natural microscopic scale of the problem, we use a 
dimensional system of units, measuring distance in Angstroms, energy and temperature in 
Electronvolts. The mass of atoms is chosen as unity: m — 1. In the remainder of the paper, 
the units of other dimensional physical quantities are omitted, but they are expressed in 
terms of the above units. 

We use exponential interactions between atoms corresponding to the repulsion between 
atomic cores, 

V mt (r) = V exp(-f3 r), (3) 

where V = 10 eV is the amplitude and inverse of /3 = 0.85 A -1 determines the range of 
interaction. 

To model the substrate, we used in the simulation the true 3D external potential, periodic 
in the (x, y) plane (with the rectangular symmetry) and parabolic in the z direction, 
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VA*; a, e, s) = ''TJZJ^ (5) 



where uo z is the frequency of normal vibration of a single atom, and 

1 (1 + sf\l - cos(27re/a)] 
2 £ 1 + s 2 - 2s cos(27rx/a) 

is the deformable Peyrard-Remoissenet potential [IS], and the parameter s, \s\ < 1, describes 
its shape. The choice of the lattice constants, a sx = 2.74 A and a sy = 4.47 A , and of the 
energy barriers, e sx = 0.46 eV and e sy = 0.76 eV, provide a high anisotropy of this potential, 
which can be viewed as the set of "channels" with the corrugated bottoms, oriented along 
the x direction. This potential is typical for the furrowed crystal surfaces, namely our 
parameters were chosen for the Na-W(112) adsystem ||10|| . The frequency u x of a single-atom 
vibration along the x direction is connected to the shape parameter s x by the relationship 
uj x = ujq(1 + s x )/(l — s x ), where ujq = (e sx /2m) 1 ' 2 (2ir / a sx ) ■ We have chosen typical values 
s x = 0.2 and s y = 0.4, which lead to the frequencies of atomic vibrations ui x = 1.65 and 
ujy = 2.02 respectively. Finally, we took u z = \{uj x + uj y ) = 1.84. Note, that our choice 
of parameters does not claim to be a quantitative interpretation of the concrete adsystem, 
because the model is oversimplified. However, we do believe on a qualitative description 
of the effect under investigation and claim that typical adsystems should exhibit similar 
behaviors. 

In the present work we study the behavior of the system in a wide range of frictions in 



the underdamped limit C w x , corresponding to typical adsystems []T6|| . In the simulation, 



we first look for the minimum-energy configuration of the system. Then, we adiabatically 
increase temperature and force and measure the mobility B for given values T and F (this 



procedure was described in detail in flID"|). In order to emphasize the phase transitions, the 
system is studied at a very low substrate temperature, T = 0.0005 eV. 

An important parameter of the generalized FK model is the atomic concentration. For 
the repulsive interatomic interaction used in the present work, we have to impose the peri- 
odic boundary conditions in the x and y directions in order to fix the concentration. Namely, 
we place N atoms into the fixed area L x x L y , where L x = M x a sx and L y = M y a sy , so that 
the dimensionless atomic concentration (the so-called coverage in surface physics) is equal 
to 9 = N/M (M = M x M y ). The atomic concentration in the FK system plays a crucial 
role since it defines the number of quasiparticle excitations, i.e. the number of geometrical 
(residual) kinks. These excitations can be defined for any background commensurate atomic 



structure 9 = p/q, where p and q are relative prime integers fll3yi7|| . If the concentration 9 



slightly deviates from the background value 9 , the ground state of the system corresponds 
to large domains with background commensurate coverage 9 , separated by localized incom- 
mensurate zones of compression (expansion) called kinks (antikinks). When the background 
commensurate coverage is trivial, 9o = 1/q, the kinks defined on this structure are called 
trivial kinks [13|]. Besides, a nontrivial background coverage #o — p/q (p / 1, with complex 
elementary cells consisting of p atoms) can be represented as a lattice of trivial kinks, defined 
on the background of the closest trivial structure. Therefore, in the latter case a deviation 
from the #o = p/q structure can be represented as a discommensuration in the lattice of 
trivial kinks, i.e. a "kink in the kink lattice", called superkink [|13 . 



As in the simulation we study finite systems, we have to chose an appropriate system 
size to insert N k kinks into the 9 = p/q commensurate background structure; the integers 
iV and M must satisfy the equation [[H]] 
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qN = pM + N k a , 



(6) 



where the topological charge a is o = +1 for the kink and a = — 1 for the antikink. As 
background structures, we discuss here two interesting cases, the trivial coverage 9 = 1/2, 
and the complex coverage 8 = 2/3. 



III. SIMULATION RESULTS. 

To simplify the problem, in the present work we only consider the quasi-one- dimensional 
case, putting M y = 1, so that all chains move in the same way (however, the interaction 
between the atoms, as well as the atomic motion still has the 3D character). Let us note 
here, that this simplified choice leads only to a minor difference in system behavior in com- 
parison with true 2D FK system with M y > 1 . Namely, the exact critical force depends 
slightly on the external conditions, which means that the transitions do not occur simulta- 
neously in all the chains. A careful examination of the behavior in different chains shows 
an enhanced transition to sliding state due to cooperative effects in the second dimension 
with an exponential law pi. On the contrary, the transition of a chain to the locked state 



is almost independent of neighboring chains. 



A. Trivial 8q = 1/2 background coverage 

First, let us consider the ground state which corresponds to domains of the trivial 6q = 
1/2 coverage separated by trivial kinks: 9 = 21/40. Namely, we took N = 105 and M x = 200, 
having thus ten kinks over the length under investigation with an average spacing of 20a sx 
between kinks. 

The generic B(F) dependence for the friction constant 77 = 0.12u x is presented in 
Fig. |l|(a). In this figure, as well as in all figures below, it will be convenient to ex- 
press force F in the units of the constant "forward" critical force for a single particle 
Ff Q = Ctte sx /ci sx w 0.607 (C ~ 1.15 for the shape of the external potential used in the 
present simulation). During the force- increasing process, one can distinguish several steps in 
the B(F) dependence, which correspond to a hierarchy of depinnings of kinks [jnj. Namely, 
at F < F t k ~ 0.23-F/o the system is in the completely immobile state, while at F > F t k the 
system has nonzero mobility B t k due to the running state of trivial kinks (Fig. a). It was 
shown in |i~0 |, that the force Ftk is related to the degradation of the periodic Peierls-Nabarro 



potential Sp^ for the trivial kinks: the kinks start to slide at F > F t k ~ K£pN/a ax . The 
second abrupt increase of the mobility to the value B m takes place when the force exceeds 
a certain threshold F pa i r w 0.35-F/o, connected with the vanishing of the energy barrier for 
creation of additional kink-antikink pairs in the system, so the number of mass carriers in 
the system increases, leading to the increase of the mobility. The remarkable tendency for 
the running kinks is to come closer to each other, thus bunching into compact groups. This 
tendency is especially enhanced after the second threshold F pair , where the concentration 
of kinks is large. The bunched kinks build up dense groups of immobile atoms with 6=1, 
while the rest of the system consists of running atoms (corresponding to the running state 
of force-excited antikinks). This state, which is very reminiscent of a traffic jam (Fig. 0(b)), 
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survives till the last threshold force Ff ~ 0.53_F/ , after which all the atoms are sliding over 
the periodic potential, and the system reaches the highest possible value of the mobility 
Bf = 1/mr). 

During the force- decreasing process, the system jumps back to the immobile state at the 
critical force Fb, and one can see a large hysteresis which, contrary to the one-particle case, 



survives at nonzero temperature |T(| 



As was mentioned in the Introduction, for a single Brownian particle the "forward" force 
FfQ is independent on the friction, while the "backward" force Fj, is linearly proportional 
to the friction. By contrast, the situation is more complex in the case of interacting atoms. 
Fig. |3| represents the dynamical phase diagram in the (F, rf) plane, where we plot the critical 
forces F t k, F pair , Ff, and Fb versus the friction coefficient rj. 

Let us consider first the forward transition from the locked to the running state. One 
can distinguish two regions of friction corresponding to different scenarios of the transition. 

(i) At very low frictions, rj < 0.05u x , there is no intermediate stages. When the force 
increases, the system jumps from the locked to the running state directly at the force Ff, 
this force being exactly equal to the critical force Ftk for the kink transition to the running 
state. 

(ii) By contrast, at larger frictions, rj > 0.05u x , all the above-mentioned intermediate 
stages with running kinks do exist. The second difference from the very low friction region 
is that the "forward" force Ff and the kink-antikink nucleation force F pair are friction- 
dependent, namely they grow with friction increasing, while the critical "kink" force F t k 
remains constant. At not too high frictions, the critical force Ff first increases approximately 
linearly with the friction 77, while at higher frictions the increase of Ff slows down, and finally 
tends to a constant value Ff, close to the value Ff for noninteracting atoms (Fig. |3|). 

Unfortunately, we failed to plot the critical force F pair in Fig. ^ at r] > 0.14w x , because 
the "traffic jam" regime is not well defined anymore at these higher frictions, probably 
because of the small number of atoms in the immobile group (Fig. 0(b)) at 9 = 21/40 for 
the finite-size system. However, at a higher concentration 9 = 21/31 (see next subsection 
MB), where the immobile groups are larger, the "traffic jam" regime is well defined at all 
studied frictions. 

The backward transition from the running to the locked state has one interesting feature. 
Namely, the "backward" critical force Fb grows with rj increasing (notice, that at not too 
high frictions, Fb{rj) exactly matches the law Fbo — 4f]y/me s /7r for noninteracting atoms, 
shown by the dashed line in Fig. |3|). If the friction is high enough, Fb is larger than the 
critical force for kinks F t k =const. Therefore, at rj < 0.1u x the system jumps back to 
the completely locked state of atoms, while at 77 > 0.1u x the backward transition has also 
multistep character: when the force decreases, it first occurs a drop of the mobility B from 
Bf to kink- mobility value B t k, and only at some F < F t k the mobility vanishes finally. 



B. Complex #0 = 2/3 background coverage 

Now let us describe the behavior of the system at the concentration 9 = 21/31 (we used 
N = 105 and M x = 155). The ground state in this case corresponds to domains of complex 
#0 = 2/3 commensurate structure, separated by superkinks with an average spacing 30a sx 
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between them. On the other hand, the 8 = 2/3 structure can be viewed as a dense lattice 
of trivial kinks defined on the background of the 9 — 1/2 structure. 

This specificity clearly manifests itself in the B(F) dependence plotted in Fig. |l|(b) for 
the same friction constant rj = 0.12u x as in Fig. |l|(a). During the force-increasing process, 
one can distinguish two sharp steps of increasing of the mobility B. The first one, at 



F ~ F sk = 0.08-F/o, corresponds to the situation where the superkinks start to slide ||10 
whereas the second step, occurring at F ~ 0.18F/ for 77 = 0. 12u x , corresponds to the 
transition of the trivial kinks to the running state, as can be seen from the phase diagram 
in Fig. |J At higher forces the increase of the mobility is smoother, and the system again 
reaches intermediate traffic jam regime with B = B m , where the mobility depends on the 
coverage as B m {6) oc Bf (1 — 9)/9, Bf = 1/mi] being the maximum mobility [|K|. 

The dynamical phase diagram (F, rj) for the complex coverage (9 = 21/31 is shown in the 
Fig. f|, where we plot F tk , Ff, F b (as in Fig. |3|) and, additionally, F sk . However, F pair is not 
plotted because the system goes to the "traffic jam" regime quite smoothly, and therefore it 
is difficult to resolve F pair . 

One can see first, that the critical force for superkinks F sk is independent on friction, while 
the force F tk increases with the friction 77 and reaches the constant value F tk ~ 0.24F/Q at 
high friction, which almost exactly matches the value F tk pa 0.23Ff for the case of 9 = 21/40 
coverage. 

The noticeable difference with the case of 9 = 21/40 coverage is in the dependence of 
the forward critical force Ff(rj). First, as the value of F sk is very low, the low- friction region 
with Ff =const is very narrow (77 < 0.01u x ), and was not resolved in the simulation. Second, 
the initial linear increase of Ff with 77 has a larger slope than for the 9 = 21/40 case; at 
high frictions the critical force Ff also tends to a constant value FJ pa Ff . Third, there 
are some irregularities in the Ff{rf) dependence at intermediate frictions. Namely, from the 
plots B(F) at higher frictions (similar to the case plotted in Fig. |l](b)) we found that in 
the "traffic jam" regime the mobility does not have a constant value B m as it is for lower 
frictions (the last plateau in Fig. [l|(b)), but at certain intervals of frictions B(F) exhibits 
an additional plateau with a larger value of the mobility, B = B m [B m < B m < Bf). This 
should result in a lower value of the critical force Ff due to the larger kinetic energy of the 
system in this intermediate state (see discussion in the next section). 

Finally, we observed that the "backward" critical force Fb is almost independent on cov- 
erage, and its dependence on friction is again well described by the corresponding expression 
for noninteracting atoms. 



IV. DISCUSSION 

The most interesting feature of both phase diagrams in Figs. [| and |] is that the criti- 
cal forces, separating various intermediate stages during the forward transition, are friction 
dependent, except the first critical force, which corresponds to the completely locked pre- 
ceding stage. From this, we may conclude that the kinetic energy of the system in the 
preceding stage defines the transition to the following stage. Below we present a discussion 
which allows to rationalize the contribution of the kinetic energy of the running kinks to the 
dynamical phase transition towards the running state of atoms. 
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First let us consider the lowest interval of frictions, 77 < 0.1oj x in the phase diagrams 
of Figs. § and [|. The explanation of the dependence of Ff on friction 77 in this region can 
be done solely with the help of kinetic arguments in the following way. Assume that there 
is a "critical" velocity of kink v c , above which the running kink cannot exist as a stable 



quasiparticle (as was found by Peyrard and Kruskal |I8 |, such a critical velocity does exist 
for the running kinks in the frictionless case). Then, if at given values of the force F and the 
friction 77 the kink drift velocity (vk) = F/mkrjk is higher than v c , the kink should destroy 
itself as soon as it starts to move, this immediately will cause an avalanche driving the whole 
system to the totally running state of atoms. 

If one makes simplifying assumption that 77^ ~ 77, the drift velocity of kinks is (vk) oc F/rj 
and, therefore, the region in the phase diagram (F,rj), where the running kinks are stable, 
is bounded by the straight line F oc v c r\. This simple linear dependence describes quite well 
the dependence Ff{r\) for 77 < Clu^ (see dotted lines in Figs. [|, f|). Thus, for instance, in 
the 9 = 21/40 case (Fig. for 77 < 0.05uj x , when the applied force exceeds the threshold F tk , 
the stationary drift velocity for the kink is higher than v c , and the system goes directly to 
the running state of atoms. This mechanism provides, that Ff = F t k = const in this region. 
However, at higher frictions, 7/ > 0.05^, the kink can move as a stable quasiparticle, and the 
transition to the running state takes place at a higher force Ff cx 7/ (Fig. ^). Such a linear 
dependence holds also (and even better) for the coverage 9 = 21/31 at not too high friction 
(see Fig. ^). Unfortunately, we failed to find the value of the critical velocity explicitly, 
because we do not know the values for the kink mass and the friction coefficient 77^. 
However, we observed that this critical velocity v c is proportional to the sound velocity c 
(defined as c 2 ~ a A V-^ t (aA) , where a a = a S x/9o is the average interatomic distance along the 
chain). Indeed, from the slopes of dotted lines in Figs. [3], [|the ratio of the critical velocities 
for both studied coverages is t>( 21 / 40 ) /t>( 21 / 31 ) « 0.72 ± 0.07, and the corresponding ratio of 
the sound velocities for the background coverages is c*- 1 / 2 - 1 / c^ 3 - 1 = 0.74. 

At higher frictions, 77 > 0.1u x , Ff starts to deviate significantly from the simple linear 
law, and finally tends to a constant value FJ. This behavior may be qualitatively explained 
if we take into account that the running kinks stimulate the transition of the system to the 
totally running state. Indeed, owing to their nonzero kinetic energy T^mk, the running kinks 
effectively reduce the average energy barrier for the transition of all the atoms to the running 
state; therefore the effective barrier e e ff should be lower than e sx . Let us approximate the 
effective barrier as e e // = £sx — Tkmk, where oc (ffe) 2 oc F 2 /rf. Then, the critical force 
for the transition to the running state, Ff = Cire e ff /a 8X , is determined by the following 
equation, 

^Ff + F f = Ff*, (7) 

where a > is a phenomeno logical coefficient, and FJ w Ff = Cne sx /a sx . Eq. (0) provides 

a qualitative agreement with the simulation data: at small 77 the forward force is proportional 

/ p* 

to 77, Ff(rj) = y -^77, while at 77 — » 00 we have Ff — ► FJ . Unfortunately, the quantitative 
agreement is not satisfactory. Even if Eq.(^) can provide a rather good fitting of the Ff(rj) 
dependence in Fig. the nonmonotonic peculiarities of Ff{rj) in Fig. |] remains unexplained. 
We recall, however, that the main simplification made above, is that the friction coefficient 
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for the running kinks rjk is equal to r\. But the running kink experiences also the intrinsic 
friction rj int , so that the correct expression is rjk = V+Vint, where r\ int depends very nonlinearly 
and resonantly on the kink velocity (v^) fT8|-pO||. The intrinsic part of friction arises because 
the propagating kink excites small-amplitude waves (phonons) in the discrete atomic chain. 
At certain drift velocities the translational motion of kink (i.e. the frequency with which 
kinks "hits" the particles) can come in a resonance with the allowed phonon frequencies, 
which results in bigger losses of energy by kink and increases effective damping. This is very 
likely why one can observe nonmonotonic variation of Ff(rj) in Fig. |j. 

The proposed qualitative picture may be applied not only to Ff, but also to other 
intermediate critical forces, such as F pair or F tk in the case of the complex background 
9 = 21/31 coverage. For example, the critical force for the trivial kinks F t k increases with 
the friction in the 9 = 21/31 case, when the preceding stage is characterized by the running 
superkinks (Fig. f|); F tk tends to the constant value found for trivial background coverage 
9 = 21/40 only at high frictions, when the kinetic energy contribution of superkinks to the 
initialization of the trivial kinks motion becomes unessential. 

Finally, let us comment the fact, that the asymptotic (n — > oo) forward critical force 
F*, for both studied coverages, has been found equal to the critical force for a single atom 
Ffo = Ctt6 sx / a sx . This is easy to understand if one takes into account that the studied atomic 



chain corresponds to a low-coupling limit [13], because the dimensionless elastic constant 

2 

g e ff = 2 ^i x Vln t (aA) is well below unity for both studied coverages, namely g e ff G [0.06, 0.2] 
for 9 G [1/2,2/3]. Therefore, the interaction between atoms should not lead to a significant 
change of the barrier for the transfer of an atom to the running state unless the running 
kinks contribute to this transition. 



V. CONCLUSIONS 

Thus, we have studied the dynamical phase transition from the locked to the running 
state of interacting atoms in a periodic external potential under the action of a dc external 
force in the underdamped limit of a generalized Frenkel-Kontorova model. This transition 
proceeds by a complex multistep scenario, which can be treated as a hierarchy of depinnings 
of quasiparticle excitations of the FK model (kinks). The interesting feature of the transi- 
tion is that the critical forces separating different intermediate stages during the "forward" 
transition, are friction dependent, except the first critical force, which corresponds to the 
transition from the completely locked ground state. This reflects the dynamical nature of 
transitions between intermediate stages, namely that the kinetic energy of running kinks in 
the preceding stage significantly contributes to the transition towards the following stage. 

On the basis of simulation results, we have proposed a phenomenological approach which 
qualitatively explains the observed friction dependences of the critical forces. According 
to this approach, for low frictions the critical quantity determining the criterion for the 
transition to the totally running state of atoms, is the drift velocity of kinks at the preceding 
stage. This approach leads to the simple phenomenological equation (|7|) for the dependence 
of the forward critical force on friction. However, for better quantitative description of the 
simulation results one has to take into account the resonant interaction of the running kinks 



with the phonons excited in the atomic chain [18— 20| . Work along this line is in progress 
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FIGURES 



FIG. 1. The mobility B versus external force F for the quasi-one-dimensional FK model: (a) 
for 9 = 21/40 coverage (trivial kinks on the background of trivial 9q = 1/2 structure); (b) for 
9 = 21/31 coverage (superkinks on the background of the complex 6q = 2/3 structure). The plus 
signs and the solid curve correspond to force increasing, the cross signs and the dashed curve to 
force decreasing. The arrows indicate the hysteresis. The force F is expressed in terms of the 
"forward" force for a single atom Ffo = CiT£ sx /a sx defined in the Introduction. 

FIG. 2. Illustration of the atomic motion in the regime of running kinks (a) and in the "traffic 
jam" regime (b). Immobile atoms are denoted with gray circles, running atoms and atoms in the 
kink regions are denoted with black circles, arrows indicate the direction of atomic motion. 

FIG. 3. Dynamical phase diagram in the (F, rj) plane for the quasi-one-dimensional FK model 
at the 9 = 21/40, i.e. for the trivial kinks on the background of simple 9q = 1/2 structure. The 
force F is scaled in the same way as in Fig. ED. 

FIG. 4. Dynamical phase diagram in the {F,rf) plane for the 9 = 21/31 coverage. The force F 
is scaled in the same way as in Fig. [p. 



12 



a) o + o + o + o + o + » •-•+ o + o + o + o + o +o 



b) 



+ 



traffic jam regime 
+•+►+ O O O O O O 



1.0 



0.8 



0.8 



0.6 



03 



F/F 



0.4 



0.2 



Running 

a) 9 2*740 t] = 0:12cj x 



atoms 




0.0 0.1 



0.2 0.3 0.4 

F/F f0 



0.5 0.6 



1.0 



0.8 



0.6 



F/F 



fO 



0.4 



0.2 



0.0 



Running 
atoms 



V9 



1^ ^ 



.4 



f 



i F 

1 



/I 



✓I 



■^z' Running 
^ kinks 




tk 

i i i i 



sk 



X X X X X X X X X X X X X X X X — X — X — X — X — X — X — X — X — X 

/z Locked state 



0.0 



0.1 



0.2 
tVco 



0.3 



0.4 



